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ABSTRACT 

Context. High resolution X-ray observations of classical T Tauri stars (CTTSs) show a soft X-ray excess due to high density plasma 
(«e = 10" - 10'^ cm"^). This emission has been attributed to shock-heated accreting material impacting onto the stellar surface. 
Aims. We investigate the observability of the shock-heated accreting material in the X-ray band as a function of the accretion stream 
properties (velocity, density, and metal abundance) in the case of plasma-/? « 1 (thermal pressure <K magnetic pressure) in the post- 
shock zone. 

Methods. We use a 1-D hydrodynamic model describing the impact of an accretion stream onto the chromosphere of a CTTS, 
including the effects of radiative cooling, gravity stratification and thermal conduction. We explore the space of relevant parameters 
and synthesize from the model results the X-ray emission in the [0.5 - 8.0] keV band and in the resonance lines of O vii (21.60 A) 
and Neix (13.45 A), taking into account the absorption from the chromosphere. 

Results. The accretion stream properties largely influence the temperature and the stand-off height of the shocked slab and its sinking 
in the chromosphere, determining the observability of the shocked plasma affected by chromospheric absorption. Our model predicts 
that X-ray observations preferentially detect emission from low density and high velocity shocked accretion streams due to the large 
absorption of dense post-shock plasma. In all the cases examined, the post-shock zone exhibits quasi-periodic oscillations due to 
thermal instabilities with periods ranging from 3 x 10"^ to 4 x 10^ s. In the case of inhomogeneous streams and/3 <g: 1, the shock 
oscillations are hardly detectable. 

Conclusions. We suggest that, if accretion streams are inhomogeneous, the selection effect introduced by the absorption on observable 
plasma components may easily explain the discrepancy between the accretion rate measured by optical and X-ray data as well as the 
different densities measured using different He-like triplets in the X-ray band. 
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1. Introduction More recently, high resolution (R ~ 600) X-ray observa- 

tions of some CTTSs (TW Hya, BP Tau, V4046 Sgr, Hen 3- 
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i^ ... ,. . , 11 J 1 ■ 1 T- T- • . 600, MP Mus and RU Lupi) have revealed a soft X-ray excess 

C_!? Young accreting stars, also called classical T Tauri stars , , , , ^ ' t, ^ -, , /^f, t^ , , 

/r-rrrrc \ u u ii . j- j j • .u i . j j i. .u produccd by plasma at temperature 7 ~ 2 - 3 X 10° K and elec- 

(CTTSs), have been well studied during the last decades, both ^ , , • ,on i,^n -^ t^ ,-i i-^„^ i 

J, ^, , . 1 . ,. .• 1 • . r • ^ .u 1. tron number density «„ ~ 10 - 10 cm ( Kastner et al.l 2002 : 

rrom the observational and theoretical point or view. On the ba- rr-, — : n UoA/i [TTt. — ; n rT/^^^ ,, , ' 7i 

. ,. , , , . . u • .• 'Schmitt et al.' 120051; iGunther et al.l '2006; Huenemo erder et all 



sis of the largely accepted magnetospheric accretion scenario 
gas from the circumstellar disk accretes onto the star surface 



2007; Areiroffi et al. 2007; Robrade & Schmitt 2007). This com- 
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V _ ^^ ^ _._... ._.... 

!_! ' fTi ''""" ""^ wiw^i.i...w,^^ ,^j^i^ ^ww^w.w., >^..... ...w ...^ ..^i^^ww ponentcouldbeproducedby the shock-heated plasma at the base 

H roUowing the magnetic held lines or the stellar magnetosphere ^„ , . ^ , , , . • , • i i • i 

Cn ] „ . if i ,r>r>i k o- .u .11 J- 1 • . . J . or the accretion column, the plasma density being much higher 

. . . 1 OKoemgll 1199 1V Since the circumstellar disk is truncated at , , , ,„in -t,-. , j- ^ „ . ^ 

^; — T^TTTT j- .• • 1 . J . 1 •.• than that (Wp < 10 cm ) measured tor coronae or active stars 

~ 3- 10 stellar radii, accretinggas IS acceleratedup to velocities r= — — — Vi u„^ .i ™ ,. j, • , ■ ■ ^ 

c .-irvrv ^rvrvi -1 i. .u .11 •. .• i£ ij j u (ITesta et al.M 2004'). The discovery or optical depths enects on 

of ~ 200 - 600 km s by the stellar gravitational held and, when ^, — ~ — r^^n^^ / , ,■ , , 

^, , . , J .^ , r.u .11 . u •. the X-ray spectrum of MP Mus further supports this hypothesis 

impacting on the high-density layers of the stellar atmosphere, it rr — ■ A ^ , innfiQh 

forms a strong shock t hat heats up the accretion colu mn to tem- — ^ 

peratures of few MK (Gullbring'l994t lLamziiilll998h . The mag- The idea that the soft X-ray excess detected in CTTSs 

netospheric accretion paradigm is supported by various observa- is due to shocks formed at the impact of accretion columns 

tional evidence: optical emission line s with inverse P Cygn i pro- onto the stellar surface has recently received a convincing the- 

files indicate infalling material (e.g. lEdwards et aP 1 19941) . UV oretical support by time-de pendent models of radiative accre- 

and optical ex cess comes from a photosphere heated by the ac - tion shocks in CTT Ss (Kol dobaet al.l 120081; ISacco etaLl 120081 ; 

cretion shock (ICalvet & Gullbriiiglll998l ; [Gullbring et alJl2000h . lOrlando et al.l 12010 ). In particular, these studies have predicted 

and transition emission lines in the far-UV are produced by post-shock plasma characterized by density and temperature 

the shock-heated material at teraperat ure of about 10^ K (e.g. values in the range observed and global shock oscillations in- 

lArdila et ani2002i iHerczeg et al.ll2005h . duced by radiative cooling instabilities analogous to those pre- 
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dieted in other astrophysical contexts (e. g. [Sutherland & Dopita 



u cont exts (e.g.laum t 
f' lQQfiHSafiedllQQa r; 



.1993; 'DoDita & Sutherland" 19961 |Safier||I998|; iSutherland et al 
,2003a b; Mignone 2005). Sacco et al. (2008; hereafter Paper I) 
developed a one-dimensional hydrodynamic model with the aim 
to investigate the dynamics and the stability of the shock-heated 
accreting material and the role of the stellar chromosphere in 
determining the position and the thickness of the shocked re- 
gion. The model takes into account all the important physical 
effects, i.e the gravity stratification, the thermal conduction, the 
radiative losses from an optically thin plasma, and a detailed de- 
scription of the stellar chromosphere. For simulations based on 
the parameters of MP Mus, they synthesized the high resolu- 
tion X-ray spectrum, as it would be observed with the Reflection 
Grating Spectrometers (RGS) on board the XMM-Newton satel- 
hte, and found an e xcellent agreement with the observations 
(lArgiroffi et al.ll2007h . 

Although the time-dependent models of radiative shocks 
support well the origin of the soft X-ray excess in CTTSs, sev- 
eral observational points remain still unclear: a) the mass accre- 
tion rates derived from X-ray observations are generally lower 
than those derived fr om optical and UV data by one or two 
orders of magnitude (ISchmitt et al.l l2005t lGiintheretaDl2007t 
lArgiroffi et al.l2009l: Curran et al. 2010, in prep.); b) some young 
accreting objects (T Tau, AB Aur, HD 163296) do not ex- 
hibit spectral features indicating dense X-ray emitting plasma 
(Gudeletal. 2007 ; Telleschi et al. 2007; Gunther & Schmitt 
l2009l) : c) a detailed time series analysis of the soft X-ray 
emission from TW Hydrae revealed no periodic variations 
dPrake et al.ll2009h . In addition, the absorption from the stellar 
atmosphere could play an important role in the observability of 
accretion sho cks in X-rays, especially for high density streams 
(lDrakell2005h . but up to now this effect has been not fully ex- 
plored. 

More recently. iBrickhouse et al.l(l2010h performed a long ex- 
posure X-ray observation of TW Hydrae, using the Chandra 
High Energy Trasmission Grating, that provides a rich set of di- 
agnostics for electron temperature, electron density and the hy- 
drogen column density. Temperature (Tg ~ 2.5 MK) and elec- 
tron density («e » 3.0 x lO'^ cm""*) derived from the He-like 
Ne IX line ratio diagnostics are in agreement with a standard ac- 
cretion shock model describing a single stream with uniform 
density and velocity. However, this model cannot explain the 
lower values of plasma density (we ~ 6.0 x 10" cm"-') and tem- 
perature (Tf. w 1.5 MK) derived from the Ovii line ratio, be- 
cause single stream models predict that the density of the post- 
shock plasma increases as the temperature decrea ses. In order 
to explain both the O vii and Ne ix diagnostics, Bric khouse et alJ 
( 1201 0) suggested that the X-ray emission originates from three 
plasma components: a hot (Tg » 10 MK) corona, a high den- 
sity (He ~ 6.0 X 10^^ cm"^ Te ~ 3.0 MK) post-shock region 
close to the shock surface, and a cold less dense («e ~ 2 x 10" 
cm~^, Je ~ 2.0 MK) post-shock cooling region, with 300 times 
more volume and 30 times more mass than that of the post shock 
region itself. 

In this paper we analyze the observability of the X-ray emis- 
sion from the post-shock plasma as a function of the properties 
of the accretion stream (density, velocity, and metal abundance). 
The main targets of our work are: a) determining the main signa- 
tures of the accretion shock in the X-ray band; b) estimating how 
the X-ray luminosity from the shock-heated plasma depends on 
the properties of the accretion stream; c) investigating the influ- 
ence of the absorption from the stellar atmosphere on the ob- 
served X-ray emission from the post-shock zone; d) understand- 
ing if the discrepancy between accretion rates measured from 



optical and X-ray data as well as the different densities mea- 
sured from the O vii and Ne ix He-like ions can be explained in 
the framework of 1-D hydrodynamic models, assuming plasma- 
/? «; 1. To answer these questions, we perform a set of 1-D hy- 
drodynamical simulations using the model introduced in Paper I. 
We explore the parameter space of the accretion stream, namely 
its mass density, velocity, and metal abundance. From the re- 
sults of the simulations, we synthesize the X-ray emission aris- 
ing from the post-shock zone in the [0.5 - 8.0] keV band and 
in the resonance lines of the He-like ions O vn and Ne ix, taking 
into account the effect of the absorption from the surrounding 
stellar chromosphere. 

The paper is organized as follow: Sect. |2] describes the nu- 
merical setup and the space of the physical parameters explored 
by the simulations; Sect. |3] reports the results focusing on the 
properties of the post-shock zone and on the X-ray emission; 
Sect. |4] discusses our results compared with the properties of 
CTTSs observed by high resolution X-ray spectrographs; in 
Sect.|5] we draw our conclusions. 

2. The Model 

2.1. The numerical setup 

We adopt the model introduced in Paper I. We assume that accre- 
tion occurs along the magnetic flux tube linking the circumstel- 
lar disk to the star and that plasma moves and transports energy 
only along the magnetic field lines. This assumption is valid for 
values of the plasma parameter y6 ^ 1 (where /3 = gas pres- 
sure / magnetic pressure) in the accretion column. We assume 
the accretion streams to be perpendicular to the stellar surface 
and focus our analysis on the portion of the stream close to the 
star with a constant mass accretion rate. This allows us to as- 
sume the plane-parallel geometry (the maximum expected ratio 
between the hot slab thickness and the stellar radius is < 0.2). A 
schematic description of the system geometry is shown in Fig. 
[1] The model is one-dimensional and describes the impact of the 
stream onto the chromosphere along the coordinate s (see Fig. 
[1). The plasma dynamics is described by solving numerically 
the time-dependent fluid equations of mass, momentum, and en- 
ergy conservation for a compressible conducting and optically 
thin plasma fluid above the stellar surface: 



dp dpu 
'dt^~ds 



0, 



dpu d(P + pu ) 



dpE d(pE + P)u 
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57 



■ nenHA(r) , 



(1) 
(2) 
(3) 



Mr - 1) 



P = [1 + a(p, T)] 






where t is the time; s is the coordinate along the magnetic field 
lines (see Fig. [U; u is the plasma velocity; p - /umunn is the 
mass density; jj is the mean atomic mass, that ranges between 
1.277 and 1.381 times hydrogen mass as function of the adopted 
metal abundance; niu is the mass of the hydrogen atom; Wg and 
riu are the electron number density and the hydrogen number 
density, respectively; P is the thermal pressure; g{s) is the grav- 
ity of a star with a mass M = 0.8 Mq and a radius R = 1.3 Rq; T 
is the plasma temperature; E - e+u /2 is the total gas energy per 
unit mass; e is the internal energy per unit mass; y = 5/3 is the 



G. G. Sacco et al.: On the observability of T Tauri accretion shocks in the X-ray band 



CO 
CD 

■4—' 

CO 

c 

O 

o 
o 



accr 
stre 

Uaoc 


3tion 
am 


corona 


T 




hot 
slab 

Ips 


chromosphere 



'abs 



Fig. 1. Schematic description of the system geometry. Ips and 
hsink are the expected thickness of the post-shock zone and the 
distance of the base of the slab from the transition region be- 
tween the chromosphere and the stellar corona, respectively (see 
Sect. 12.21 ). Sahs is the threshold below which, we assume the X- 
ray emission is fully absorbed (see Sect. 12.4b . 



ratio between specific heats; k^ is the Boltzmann constant; q is 
the heat flux; A(r) is the radiative losses function per unit emis- 
sion measure; a{p, T) is the fractional ionization (rie/nu), that 
has been derived from a modified Saha equation for the solar 
chromosphere conditions (lBrown|[l973l) ; En is a parametrized 
chromospheri c heating function (En = for T > 8 x 10^ K) 
defined, as in iPeres et aP (Il982l) . to keep the unperturbed chro- 
mosphere in stable equilibrium. 

The thermal conduction includes a smooth transition be- 
tween the classical and the saturated conduction regime, as in 
lDalton&Baibusl(fT99l 



The equations a re solved numerically using the flash code 
(iFrvxell e t al.' '2000^ with its implementation of the PPM al- 
gorithm (Colella & Woodward 1984). flash is a multi-physics 
code for solving astrophysical problem s that handles an adap - 
tive mesh by the paramesh algorithm dMacNeice et al.l 120001) . 
The code has been extended with an additional module for the 
evolution of the fractional hydrogen ionization (ng/nn)- 

The extension of the computational domain depends on the 
specific parameters of the simulations (see Sect. 12.2b . The largest 
domain considered here extends over a range D = 2.4 x 10'" 
cm above the stellar surface, for the simulation that analyzes an 
accretion stream density «e = lO"' cm""*, velocity u^^c = 600 
km s"', and solar abundances (see Sect. 12.21 for a detailed de- 
scription of the space of the physical parameters explored). 

The spatial resolution adopted for each simulation depends 
on the thickness of the post-shock zone that spans a range of 6 
orders of magnitude (See Sect. 13.1b . At the highest resolution 
(stream with density Wg - 10'"' cm"-', velocity Mace = 200 km s"' 
and metal abundance ^ = 5, in solar units) we allow for a maxi- 
mum of 16 levels of refinement in the paramesh algorithm, with 
resolution increasing twice at each refinement level, and with the 
refinement criterion following the changes in density and tem- 
perature. This grid configuration yields an effective maximum 
resolution of » lO' cm at the finest level. For the simulation re- 
quiring the highest spatial resolution, we analyzed the eff'ect of 
resolution on the model solution by considering additional sim- 
ulations which use an identical setup but with a larger number 
of refinement levels. We checked that the resolution adopted in 
this paper is the best compromise between accuracy and compu- 
tational cost and that the system evolution is well described in 
its details. 

As initial condition we consider a stream with uniform den- 
sity, temperature and velocity impacting onto a hydrostatic chro- 
mosphere (see Fig. |2]l. The temperature of the stream is the 
same in all the simulations performed and does not influence 
the system evolution, because the thermal energy of the stream 
is negligible with respect to its kinetic energy. The chromo- 
sp here of the young star is described adopting the solar models 
in lVernazza et al. ( 1973), rescaled to have, at its base, a pressure 
(w 9 X 10"^ dyn cm"^) larger than the ram pressure of the streams 
considered here and an extension of 3 x lO" cm. As boundary 
conditions, the values of density, temperature, and velocity are 
fixed to the initial ones both at the top (i.e. a constant accretion 
rate) and at the base of the computational domain. We checked 
that the adopted chromosphere and boundary conditions prevent 
boundary eff'ects on the system evolution. 



1 1 

q^\ — + — 



(4) 



where ^spj - -K{T)dT jds is the classical conductive flux 
(ISpitzeil ll962) (with k(T) = 9.2 x lO^^r^/^ e rg s"' K'' cm'') 
Q'sat = -sign jdT l ds)^pcl is the satur ated flux (ICowie & McKed 
Il977l) . < 1 (Bor kowski et aljfl 9891 and references therein) and 
Cs is the isothermal sound speed. 

The radiative losses per unit of emission measure MJ) 
have been calculated w ith the PINTofALE spectral code 
(iKashvap & Drakd l200(A and the APED VI. 3 atomic line 
database (Smith etal. 2001). Since the radiative losses in the 
temperature range lO*" < T < 10^ K are dominated by the emis- 
sion lines from heavy ions, the radiative cooling depends on the 
metal abundance, (, which is one of the parameters explored in 
this work. 



2.2. Space of the physical parameters 

We explore the space defined by three parameters: the veloc- 
ity and density of the pre-shock accretion flow, and the metal 
abundance of both the stellar and the accreting matter, that are 
assumed to be equal. Note however that a discrepancy between 
chromospheric and accretion metal abundances would not af- 
fect the results of our simulations, because, as discussed in 
ISacco et aP (l2008h and in the Sect. l3.1l of this paper, the chro- 
mosphere acts as reservoir of matter stopping the accretion flow, 
but it does not influence the properties of the post-shock zone. 
Instead, diff'erent chromospheric abundances can slightly affect 
the absorption from the chromosphere as discussed in Sect. 12.41 
We consider three pre-shock velocities (m^cc - 200, 400, 600 
km s"'), 4 pre-shock densities (nacc = lO'", lO", lO'^ lO'^ 
cm"^) and 3 metal abundances (4" - 0.2 , 1.0, 5.0 in solar units re- 



siT 



trieved from i Anders & Grevessd 19891) . Both densities and metal 
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Fig. 2. Density (dashed line) and temperature (dotted line) pro- 
files for the initial unperturbed stellar chromosphere {s < 3.05 x 
10** cm) and an accretion stream (s > 3.05 x 10** cm) with hydro- 
gen number density Wacc = 10'^ cm^^ and solar metal abundance. 



Table 1. Range of relevant physical parameters explored. 



Parameter 


Range 


Units 


"ace 


10'" - 10" 


cm-^ 


i<acc 


200 - 600 


km s"' 


( 


0.2-5.0 




7-p. 


5.8 X 10" - 5.6 X 10^ 


K 


1"cool 


8.3 X 10-3 - 1000 


s 


/ps 


4.2 X 10^-1.5x10'° 


cm 


p 

^ ram 


35- 7.8 X lO-* 


dyn cm"^ 


fhink 


1.3 X 10** - 3.5 X 10"* 


cm 


ySps 


8.6 X lO-'' - 2.0 





n^cc is the stream hydrogen number density; u.^^^. is the stream velocity; 
^ is the heavy element abundance; Tps is the post-shock temperature; 
Tcooi is the radiative cooling time; /ps is the expected thickness of the 



post-shock zone; P^^ 



is the ram pressure of the accretion 



flow; /jsini; is the sinking of the post-shock zone in the chromosphere; 
ySps = ^ram/(fi^/8;r) is the plasma parameter in the post-shock zone, 
assuming a magnetic field B = I kG. 



abundances span larger ranges than those obtained from high 
spectral resolution X-ray observations. We focused our analysis 
on a large space of parameters for two main reasons. Probably 
the most important is that current high spectral resolution X- 
ray observations of CTTSs are limited to a very small sample 
of very close and bright CTTSs. These stars are older and ac- 
crete to a lower rate than typical CTTSs and it is most likely that 
physical and chemical properties of the cor responding accre- 
tion stream may differ too. Indeed. lBarv et al.l (12008) and lMartinI 
(1199 6) estimated pre-shock densities Wacc = 10'^ - 10' 
while pre-shock densities derived from X-ray data are n^ 



^'"^ cm ^, 



10' 



10 cm . Moreover, densities and metal abundances 



are measured from the spatially integrated emission including 
both the coronal and the accretion component. Typical coro- 
nal densities («e < 10'" cm^^) are much lower than densi- 
ties expected in the accretion stream and metal abundances are 
lower than photosph eric abundances observed in young stars 
(iTelleschi et al.ll2007h . Therefore, both accretion stream densi- 
ties and abundances derived from X-ray observations could be 
underestimated. 

A h euristic model that assumes the strong shock approxi- 
mation (IZel'Dovich & Raizerlll967h . stationary conditions and 
radiative cooling has been previously used to describe accre- 
tion shock physics in CTTSs (Lamzin 1998; Calvet & Gullbring 
ri998l) and to derive mass accretion rates from the X-ray emis- 
sion (lArgiroffi et al.lf2007h . This model provides post-shock re- 
gion characteristics: 



^ps 



'^ps — ^^acc : 



^cool 



1 p 

y - 1 nenHA(r) 

'ps 

Across ~ = ^cool ■> 

"ns 



P'^ ~ T,2 k '^'^ l.'t X lU Ma„ 

11 J. I-I,,J-,-, 

'ps = ^cool^ps — A-/ X lU — 

^ "ace 



1 T^l^ 
2.5xl0^i P^ 



^ "ps 



(5) 
(6) 
(7) 
(8) 
(9) 



where Mps and Wps are the post-shock velocity and density. Mace 
and Hacc are the stream pre-shock velocity and density, Tcross 
is the crossing time of the accreting material through the post- 
shock zone (which is expected to be equal to the radiative cool- 
ing time Tcooi), 7"ps is the post-shock temperature, and /p, is the 
thickness of the post-shock zone, and where the radiative cool- 
ing function has been_a£2roximatedas A(r) «i 1.6xl0"'^^r"'^^ 
erg s ' cm^ (e.g. lOrlando et alj|2005h . 

Table [T]reports the ranges of values of all the relevant phys- 
ical parameters of the simulations explored here. The first three 
rows report the three independent parameters of our analysis (i.e. 
density, velocity, and metal abundance of the stream), whereas 
the following rows report the parameters derived from the in- 
dependent ones, namely the post-shock temperature Tp, (Eq.[8ll, 
the radiative cooling time Tdoss (Eq. |6]l, the slab thickness Zp, 
(Eq. |9] and Fig. [T]i, the ram pressure of the accretion stream 
Aam = PaccMacc the sinking of the slab in the chromosphere 
/isink (i-C- the distance of the base of the slab from the transi- 
tion region between the chromosphere and the stellar corona, 
see Fig. [U, and the plasma parameter of the post-shock zone 
P - Pram/(B^/^^), whcrc the magnetic field strength is as- 
sumed to be B ~ 1 kG (which is the order of magnitude of the 
field s trengths derived from observations of CTTSs; Johns-Krulll 
l2007h . Note that the plasma /? is «: 1 for most of the accretion 
streams considered in this work (thus justifying the low-j8 as- 
sumption of our model) except for 5 cases, where B ~ 1 -2 (for 
the magnetic field assumed). On the other hand. lOrlando et aP 
(.2010.) showed that, for /3 ~ 1-5, the evolution of radiative 
shocks in 2D MHD models (thus including an explicit descrip- 
tion of the ambient magnetic field) is analogous to that described 
by 1 -D hydrodynamic models, although the amplitude of oscilla- 
tions is smaller and the frequency higher than those predicted by 
1-D models. For the 30 different cases analyzed here, the dura- 
tion of the simulations ranges from 100 to 19000 s. The duration 
of each simulation was set much longer than the initial transient 
effect, at the impact of the stream onto the chromosphere, and to 
include a large number (at least three) of oscillation periods of 
the post-shock zone (see Paper I). 
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2.3. Synthesis of the X-ray emission 

For each simulation, we synthesize the X-ray luminosity (Lx) in 
the [0.5 - 8.0] keV energy band and in the resonance lines of two 
He-like ions, namely the Ovii at 21.60 A (Lovii) and the Neix 
at 13.45 A (LNeix), which can be measured with the high reso- 
lution spectrographs on board the XMM-Newton and Chandra 
satellites. In particular: 

1. From the spatial distributions of density and temperature 
(output of the numerical simulations), we compute the emis- 
sion measure of the plasma in the /th domain cell as emj - 
n^.AstrAij, where nej is the particle number density in the 
cell, Asj is the length of the jth cell along the stream, and 
Asti is the cross section of the stream. In this paper we as- 
sume Astr = 5 X 10^" cm^, (i.e. ~ 0.5% of the stellar surface); 
for the values of stream velocity and density explored here, 
this cross section produces mass accretion rates spanning the 
range of values derived from X-ray observations of CTTSs 
(see Sect. 13.31 1. We then derive the distribution of emission 
measure versus temperature, EM(r), by binning the emis- 
sion measure values into slots of temperature; the range of 
temperature [4 < logr(K) < 8] is divided into 81 bins, all 
equal on logarithmic scale (Alog(r) = 0.05). 

2. From the EM(r) distributions, we synthesize the X-ray 
emission in the [0.5 - 8.0] keV band and in the resonance 
lines of Ovii and Neix, considering the appropriate metal 
abundance ^, using the PlNTofALE spectral code with the 
APED VI. 3 atomic line database. 

3. Finally, we derive time-average luminosities Lx, iovii, i^Neix 
over a selected time interval including, at least, three oscilla- 
tion cycles and outside initial transients. 

2.4. Absorption from the stellar atmosphere 

As pointed out by iDrakd (l2005h . a large fraction of the X-ray 
emission due to the accretion shock may be absorbed by the gas 
present in the surrounding stellar atmosphere or by the accre- 
tion column itself. The absorption from the stellar atmosphere 
mainly depends on the location of the emitting plasma. For a 
given shocked slab of material rooted in the chromosphere, pho- 
tons emitted from plasma located at the base of the slab (deep 
into the chromosphere) travel through higher density gas layers 
than photons from plasma located in the shallower portion of the 
slab. The deepness of the emitting plasma in the chromosphere 
depends on the thickness of the slab /ps and on the sinking of 
the slab in the chromosphere /isink (Sect. 12.2b to the position at 
which the ram pressure of the post-shock plasma equals the ther- 
mal pressure of the chromosphere. As discussed in Sect. 12.21 our 
simulations explore a wide range of values for /ps and /zsink (see 
Table[B. 

In order to estimate the effect of the absorption on the ob- 
served X-ray emission, we calculate the X-ray luminosities by 
assuming that only the emission from plasma located in th e shal- 
lower layers of the shocked slab can be observed (see also lDraka 
12005 ' for a previous application of this method to estimate ab- 
sorption effects). Specifically, we assume that all the emission 
due to plasma located at s < Sabs is fully absorbed, while the 
emission produced from plasma located at s > iabs is fully trans- 
mitted. The threshold ^abs is defined as the height, in the unper- 
turbed stellar chromosphere, at which the overlying atmosphere 
absorbs 50% of the energy of a trial X-ray spectrum (see also 
Fig. [TJ. We considered as trial X-ray spectrum that produced by 
a plasma at T = 1 MK. The plasma temperature does not in- 



Table 2. Values of the threshold ^abs as a function of the energy 
band and of the metal abundance. 



i 



[0.5 - 8.0] keV 
(10** cm) 



Ovii 

(10* cm) 



Neix 
(10** cm) 



0.2 
1.0 
5.0 



3.56 
3.72 
3.81 



3.57 
3.74 
3.83 



3.17 
3.43 
3.58 



fluence the threshold .Sabs used for calculating the Ne ix and the 
O VIII resonance line luminosities, while it slightly influences the 
thresholds used for the [0.5 - 8.0] keV band. However, we found 
that the results of only a few of the parameter configurations 
considered here, specifically those leading to Ips ~ hsi„k, slightly 
depend on the particular choice of the trial X-ray spectrum. For 
these configurations a more detailed description of the absorp- 
tion effect, including the dependence on the wavelength, a more 
smooth transition between the full absorption and the full trans- 
mission and a dependence of the angle of view should be per- 
formed to give a more precise answer to the observability issue. 

The values of ^abs for the three X-ray bands and for the three 
metal abundances are reported in Table |2] The value of .Sabs ob- 
viously depends on the chromospheric metal abundances, that 
we assume to be equal to the abundances of the accreting mate- 
rial. However, chromospheric abundances may be different from 
accretion stream abundances. For instance, if the chromospheric 
abundance is higher or lower than the accretion stream abun- 
dance, then we underestimate or overestimate the absorption ef- 
fect, respectively. A better knowledge of metal abundances of 
the different components of the star-disk system is required to 
better address this issue. 

We found that iabs does not change significantly if we con- 
sider the absorbed X-ray emission in the [0.5 - 8.0] keV band or 
in the O vii resonance line, because the absorption in the softer 
part of the [0.5 - 8.0] keV band (which is the same of the O vii 
line) is dominant. On the other hand, we found a significantly 
lower value of iabs in the Ne ix resonance line, this line form- 
ing at higher energy (where the absorption is lower) than O vii. 
Therefore, we consider one set of thresholds ,Sabs for the synthe- 
sis of the luminosities in the [0.5 - 8.0] keV band and in the O vii 
resonance line, and another set of ^abs for the luminosities in the 
Ne IX resonance line (see Table|2]i. 



3. Results 

3.1. Physical properties of the post-shock zone 

In all the simulations, during the first 100 - 200 s, the system 
follows the same evolution as described in Paper 1. Figure |3] 
shows, as an example, the evolution of plasma temperature, den- 
sity, pressure and velocity for a stream with «acc - 10" cm"^. 
Mace = 400 km s"' and ( - 1.0. The accretion stream penetrates 
the chromosphere and generates a transmitted shock and a re- 
verse shock. After this transient phase, the chromosphere stops 
the flow penetration where the ram pressure of the post-shock 
plasma equals the thermal pressure of the chromosphere. The 
reverse shock progressively builds up a nearly isothermal slab at 
the post-shock temperature (heating phase; left panels of Fig.[3]l. 
During this phase the intensity of the radiative cooling at the base 
of the slab gradually increases. The heating phase ends when the 
radiative cooling triggers a thermal instability that robs the post- 
shock plasma of pressure support, causing the material above 
the cooled layer to collapse back (cooling phase; right panels 
of Fig. 13). Consequently the reverse shock moves downwards to 
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Fig. 3. Time evolution of plasma temperature 
(A), density (B), pressure (C) and velocity (D) 
during a complete oscillation of the post-shock 
zone in the case of an accretion stream with 
density «acc = 10" cm"^, velocity Wacc = 400 
km s"' and metal abundance f = 1.0. The left 
(right) panels show the time evolution during 
the heating (cooling) phase. The dashed lines 
describe the profiles of the unperturbed initial 
stellar atmosphere (chromosphere and corona). 
The shaded areas mark the regions below the 
thresholds Jabs used to estimate the effects of 
absorption on the luminosity in the O vii (pale 
grey) and Neix resonance lines (see Sect. 12.4b . 



the chromosphere reducing the post-shock zone thickness. After 
the hot slab has disappeared, a new slab is re-built by the re- 
verse shock, starting a new cycle of quasi-periodic shock oscil- 
lations. The post-shock zone is separated from the cold chro- 
mosphere by a very steep transition region located where the 
chromospheric thermal pressure equals the ram pressure of the 
post-shock plasma. The thermal conduction drains energy from 
the shock-heated plasma to the chromosphere through this tran- 
sition region, acting as an additional cooling mechanism (see 
also Orlando et al. 2010). The effects of the heat conduction are 
the largest in the simulations with the highest post-shock tem- 
peratures (i.e. the highest stream velocities). 

Figure |4] and Table [3] show the main physical parameters 
characterizing the structure of the post-shock zone after the ini- 
tial transient phase, namely the maximum extension of the post- 
shock slab (Imax), the oscillation period (Pose), and the emission- 
measure-weighted temperature (Tps) of the post-shock zone. In 
order to define /max, we consider all the material hotter than 
T - 3 X 10^ K as shocked plasma. The steepness of the transition 
region makes l,„a_x poorly sensitive to this choice. The periods of 
the oscillations have been determined through a Fourier analysis 
of the length of the post-shock zone as a function of time. The 
temperature Tps has been calculated from the emission measure 



distribution EM(r) (see Sect. 12.31) . considering only bins with 
logr(K)>5. 

The maximum extension of the post-shock zone ranges be- 
tween a factor 0.3 and 1.7 of the length estimated using Eq. |9] 
whereas the post-shock temperature ranges between a factor 0.7 
and 1 .0 of the temperature estimated using Eq. [8] The simulated 
post-shock temperatures are slightly lower than those derived 
from Eq.[8]because the former are averaged over the whole post- 
shock zone, whereas the latter are estimated at the shock front, 
at the maximum extension of the slab. In addition, the simulated 
post-shock temperatures are averaged over several shock oscil- 
lations, including both the heating and the cooling phase of the 
evolution. During the former, the temperature is slightly higher 
than that derived from Eq. [8] because the plasma velocity in the 
reference frame of the shock is higher than u-^^c', the opposite is 
true during the cooling phase. 

The post-shock temperature depends on the accretion stream 
velocity, whereas both the oscillation period Pose and the length 
of the post-shock zone /max depend on the whole set of model pa- 
rameters (velocity, density and metal abundance of the stream). 
Figure |4] and Table [3] show that Pose and /^ax both span 5-6 or- 
ders of magnitude for the range of parameters explored here. In 
the complex magnetospheric accretion scenario, the mass den- 
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Table 3. Physical properties of the post-shock zone and X-ray luminosities 
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the 0.5-8.0 keV band not considering the absorption. 

the 0.5-8.0 keV band considering the absorption from the stellar atmosphere as discussed in the text. 

the O VII resonance line at 21.6 A not considering the absorption. 

the O VII resonance line at 21.6 A considering the absorption from the stellar atmosphere as discussed in the text. 

the Neix resonance line at 13.45 A not considering the absorption. 

the Ne ix resonance line at 13.45 A considering the absorption from the stellar atmosphere as discussed in the text. 



sity and velocity of the flow could v ary along the accretio n 
stream cross section (see for instance iRomanova et al.l 120041) . 
Consequently, in the case of/? « 1, an accretion stream can- 
not be described by a single 1-D model but has to be con- 
sidered as a bundle of independent fibrils, each described in 
terms of a different 1-D model, and each independent on the 
others (with different instability periods and random phases of 
the shock oscillations) due to the strong magnetic field which 
prevents mass and energy exchange across magnetic field lines. 
We expect therefore that accretion streams consisting of many 
(200 - 300) different fibrils with different instability periods and 
random phases would produce no periodic variations of the X- 
ray emissio n from shocked pla sma, as recently found in the case 
of TW Hya (iDrake et al.ll20()9l) . 

As discussed in Sect. 12.41 the key parameters for understand- 
ing the effects of the absorption by the stellar chromosphere on 
the X-ray emission from shocked accreted plasma are the thick- 
ness of the shocked slab /^ax and the sinking of the slab in the 
chromosphere to the position /isink at which the ram pressure of 
the post-shock plasma equals the thermal pressure of the chro- 
mosphere. Figure |4] shows that /max is anti-correlated with the 
density of the accretion stream and with the metal abundance, 
and is correlated with the stream velocity (see also Eq. |9]l. On 



the other hand, /isink is correlated with both the stream density 
and velocity, the ram pressure of the post-shock plasma being 
^ram '^ PaccMacc ^^^ '^^c thermal prcssurc of the chromosphere 
decreasing with height. Consequently, the shocked slab from 
high density streams with high metal abundances is expected 
to hardly emerge from the dense chromospheric layers, where 
the absorption is strong, because the post-shock zone is rooted 
deeply in the chromosphere (/isink is large) and the slab is rather 
thin (/max is small). On the other hand, streams with high veloc- 
ity form extended post-shock zones (being /^ax ^ "ace) ^^^ the 
shocked slab may easily emerge above the chromosphere even if 
it is deeply rooted in the chromosphere, the ram pressure being 
proportional to the square of velocity. 



3.2. Distribution of emission measure versus temperature 

The distribution of emission measure versus temperature EM(r) 
of the shock-heated plasma is a useful source of information 
of the plasma components contributing to the X-ray emission 
and is directly comparable with EM(r) distributions derived 
from X-ray observations (see, for instance, Argiroffi et al. 200^. 
Figure|5]shows how the time-averaged EM(r) of the post-shock 
plasma varies as a function of velocity (upper panel), density 
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Fig. 4. Periods of the shock oscillations (Poac, upper panel), max- 
imum extension of the post-shock zone (/max, middle panel) and 
the emission-measure-weighted temperature (Tps, bottom panel) 
of all the cases investigated. The labels on the x-axis of the bot- 
tom panel indicate the densities of accretion stream. Different 
symbols indicate different velocities: Mace = 200 (circle), 400 
(squares), and 600 (triangles) km s"'. Different colors indicate 
different metal abundances: ^ - 0.2 (black), 1.0 (red) and 5.0 
(green). Note that data points have been spread over each den- 
sity value to prevent overlapping of the symbols. 



(middle panel), and metal abundance (lower panel) of the ac- 
cretion stream. The case with density «acc - 10'^ cm"-', veloc- 
ity Mace = 400 km s"' and solar abundance ^ - I (red dotted 
line in Fig. |5]l is the reference in all the panels. In all cases, 
our model predicts a monolithic distribution of emission mea- 
sure of the post-shock plasma that covers the entire range of 
values below the temperature of the shock front. The ascend- 
ing part of the EM(r) distribution corresponds to cooled post- 
shock plasma of the slab and to the transition region between the 
shocked chromosphere and the hot slab. This characteristic of 
the EM (r) distribution cannot be reproduced by heuristic mod- 
els (e.g. lLamzin|[l998HCalvet & Gullbring|[T99i lArgiroffi et all 



I2007[) and has been prove d to be in agreement with observations 
re.g. lArgiroffietal.ll2009h . 
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Fig. 5. Time-averaged emission measure distribution of the post- 
shock plasma for accretion streams with the same metal abun- 



dance (^ = 1.0) and density («ac 



10'^ cm""*) but different 



velocities (upper panel), for streams with the same metal abun- 
dance (^ = 1.0) and velocity (Mace = 400 km s"') but different 
densities (middle panel), and for streams with the same density 
("ace - 10'^ cm"-') and velocity (Mace - 400 km s"') but different 
metal abundances (lower panel). 



The total emission measure of plasma above 1 MK and the 
maximum temperature of the emission measure distribution de- 
pend on the velocity of the accretion stream (upper panel in 
Fig. |5j. In fact, the temperature of the post-shock plasma de- 
pends on the square of the stream velocity (see Eq. |8]l, lead- 
ing to a shift of the EM(r) profile towards higher temperature 
for higher value of Mace- The thickness of the post-shock zone 
strongly increases with the stream velocity (see Fig.|4]and Eq.|9]l, 
so that its volume and therefore the overall emission measure of 
the shock-heated plasma increases with Mace- 

Not surprisingly, the density of the stream determines the 
overall emission measure in the hot slab (see middle panel in 
Fig. |5]i: the higher the stream density, the higher the emission 
measure of the slab. Note however that, since /ace °^ "ace (^66 Eq. 
|9]l, the emission measure of the post-shock zone depends only 
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linearly on n^cc (not on its square). Figure |5] also shows that the 
ascending branch of the EM(r) distribution is steeper for larger 
values of density. 

Finally, the EM(T) distribution depends on the metal abun- 
dance: the higher the value of ^, the lower the emission measure 
of the slab. In fact, the length of the post-shock zone depends 
inversely on the metal abundance (see Eq. |9l) as a consequence 
of more efficient radiative cooling for larger (. Consequently, the 
emission measure of the hot slab decreases for increasing (. 

3.3. X-ray emission 

X-ray luminosities derived for each simulation are reported in 
Table|3]and in the Fig.|6] where they are plotted against the mass 
accretion rate derived assuming the same accretion stream cross 
section (Asu- = 5 x 10^" cm-) used for deriving emission measure 
distribution (see Sect. 12.3b . The left panels show the results ob- 
tained not considering absorption due to the stellar atmosphere, 
whereas right panels show the results including the absorption 
effects using the method discussed in Sect. 12.41 

The dependence of the X-ray luminosities on the model 
parameters can be easily explained by interpreting the results 
shown in the right panels of Fig. |6] in the light of the EM(r) 
distributions shown in Fig. |5] and of the relations between the 
main properties of the post-shock zone and the stream param- 
eters in Fig. m The spread in the X-ray luminosity due to the 
metal abundance is less than 0.5 dex even with the large range (a 
factor 25) of metal abundances considered in this work. This re- 
sult is due to two different effects working in opposite directions. 
The EM(T) of the post-shock plasma is anti-correlated with the 
metal abundance (bottom panel of Fig. |5]l, but higher metal abun- 
dances trigger higher emission in the soft X-ray band, because 
soft X-ray emission is mainly due to line emission produced by 
heavy ions. 

The X-ray luminosity is strongly correlated with the accre- 
tion stream velocity: X-ray luminosities derived for streams at 
velocities u^^c - 400 km s"' and Mace = 600 km s"' differ 
for about one order of magnitude, whereas X-ray luminosities 
produced by streams with Mace = 200 km s"' are more than 
three orders of magnitude lower than the others. In addition, the 
spread between the streams at 400 and 600 km s"' is much less 
in the O vii resonance line (left-central panel of Fig. |6]l than in 
the [0.5 - 8.0] keV band. Evidently, the relation between X-ray 
luminosity and stream velocity strongly depends on the profile 
of the EM(r) distribution. In fact, low velocity (200 km s"') ac- 
cretion streams cannot heat up plasma to temperatures greater 
than 1 MK (top panel of Fig. |5]l, so that the accretion energy is 
emitted mainly in the extreme ultraviolet band and very faint X- 
ray emission is expected. The peak of the EM(r) for the stream 
at 400 km s"' (log T(K) - 6.4) is very close to the peak of the 
emissivity function of the O vii resonance line and it is compara- 
ble to the EM value at log r(K) = 6.4 for the case at 600 km s"' , 
explaining the small difference between the Ovii luminosities 
for streams at 400 and 600 km s '. 

The comparison between the X-ray luminosities reported in 
the right and left panels of Fig.|6] as well as the relation between 
the ratio Imax/hsmk and the model parameters (Fig|7]i, allow us to 
highlight the role of absorption from the stellar atmosphere on 
the observability of the post-shock zone in the X-ray band as a 
function of the accretion stream properties. Specifically, we find 
that: 

- Low density (wacc ^ lO" cm"-') and high velocity (Mace ^ 
400 km s"' ) streams form shocked slabs with X-ray emission 



10''^ cm ■') streams are less evident. 



es- 



poorly absorbed by the surrounding stellar chromosphere, 
the thickness of the post-shock zone being much larger than 
the sinking of the stream in the chromosphere (/^ax 3> /isink), 
so that most of the emitting plasma is located in the shal- 
lower and low density portion of the chromosphere. 

- High density («aee ^ lO'^ cm"^) and low velocity (Maee ^ 200 
km s"') accretion streams form shocked slabs with X-ray 
emission strongly absorbed by the chromosphere, being in 
these cases /ps «; /isink- 

- The effects of the absorption in the case of intermediate den- 
sity («aee ~ lO" 

pecially for those cases with /^ax ~ /^sink 
gest that except for the stream with Maee = 600 km s 
and metal abundance ^ = 0.2, all the streams with density 
"ace ~ 10'^ cm"-' generate shocked slab whose X-ray emis- 
sion is strongly absorbed (see Fig.|7ll. However, these inter- 
mediate cases probably require a more accurate analysis of 
the absorption effect that considers: a) the absorption depen- 
dence on the wavelength; b) a detailed model describing the 
absorption from the stellar atmosphere with a more smooth 
transition between the optically thin and the optically thick 
case and different angles of view; c) different metal abun- 
dances in the accretion stream and in the surrounding atmo- 
sphere. 



Our results sug- 
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Fig. 7. The ratio Imaxlhsink as function of the mass accretion rate. 
Filled and empty symbols indicate if X-ray emission from the 
stream is observable or absorbed, respectively. Symbols and col- 
ors as in Fig. |4] 



It is worth to emphasize that if accretion streams are not 
uniform in density, the chromospheric absorption triggers a se- 
lection effect, absorbing preferentially the X-ray emission from 
high density plasma components. In addition, since the chromo- 
spheric absorption does not equally affect the different density 
tracers (e.g. Ovii and Neix He-like triplets; see right panels in 
Fig.|6l), we expect that, in general, the results of the use of these 
tracers may differ significantly (see also discussion at the end of 
Sect.©. 



4. Discussion 

Observational data and 3D magnetospheric models of the star- 
disk interacti£nhi_j^oung stars draw a very complex scenario 
(see lBouvier et alJ2007l for a review of the recent results), where 



G. G. Sacco et al.: On the observability of T Tauri accretion shocks in the X-ray band 



1) 
en 
CD 






CO 
0) 



,32 



i30 



-,28 



^30 



-,26 



k 

■ 1 
ff 

t 

f 

t 

k 

■ 1 . 



^Q-11 ^Q-10 ^Q-9 ^Q- 



^Q-11 ^Q-10 ^Q-9 ^Q- 



Mccc (Mq yr") 



^ccc (Mo yr") 



Fig. 6. Luminosities in the [0.5 - 8.0] keV band 
(top panels) and in the resonance lines of the 
Ovii (middle-panels) and Neix (bottom pan- 
els) as function of the mass accretion rates. 
Absorption from the stellar chromosphere is 
considered only for the results reported on right 
panels. Symbols and colors as in Fig.|4] 



circumstellar gas accretes onto the star through many streams 
observed from different angles of view and, likely, non-uniform 
in density. Our modeling focuses on the portion of the stream 
impacting on the chromosphere where a shock develops and as- 
sumes the stream with a constant mass accretion rate. It is worth 
to emphasize, however that the model results presented here can 
be combined together in a tool-box like fashion, if the plasma- 
y6 «; 1 ( assumption which is valid for many CTTSs as discussed 
in Sect. 12.2b . allowing us to build up inhomogeneous streams 
as bundles of independent fibrils, each describable in terms of a 
different 1-D model (see discussion in Sect. 13.1b . This method 
can be easily applied to investigate the X-ray emission expected 
from a 3D density structured accretion flow, adding together the 
contribution to the emission of the distribution of fibrils each 
with its own density. In the following, we apply this idea to very 
simplified cases. In the future the method could be applied to 
describe more complex configurations. 



The hypothesis that the soft X-ray emission from CTTSs is 
mainly due to post-shock accreting plasma is supported by the 
evidence that the emission originates from high density plasma 
(we > 10" cm""*). This evidence has been found in all (six) 
CTTSs with mass lower than 1.2 Mq observed to date by high 
resolution X-ray spectrographs. In addition, no signatures of 
high density plasma have been revealed in current X-ray ob- 



servations of weak- line T Tauri stars (e.g. iKastner et al.ll2004t 
lArgiroffi et al.ll20()5h in which accretion is known to be absent. 

However, although the soft X-ray excess in CTTSs can be in- 
terpreted as due to accretion shocks, some observational results 
are still not easily explained in terms of post-shock accreting 
plasma: a) if we assume to observe all the energy emitted by 
the post-shock zone, the accretion rates derived from X-ray data 
are generally about 1-2 or ders of magni t ude lower than those 
derive d from optical data (ISchmitt et al.l 120051: iGiinther et al.l 
2007] lArgiroffi et ai] l2009t Curran et al. 2010, in prep.); b) 
plasma densities derived from the Ovii triplet a re generally 
lower than those derived from the Ne ix triplet (Brickhouseet al.l 
2010), although we expect the opposite result as the Neix 
triplet forms at higher temperatures than O vn where accretio n 
shock models predict lower densities (e.g. ISacco et al.l [20081) : 
c) no signatures of high density plasma have been found in the 
high mass (M > 2M (7)) accreting objects, namely T Tau, AB 
Aur and HP 163296 (iGiid el et al.ll2007l iTelleschi et aD l2007l 
IGiinther & Schmitlll2009l) ! 

The above issues can be addressed in the light of the results 
discussed in Sect. [3] The comparison of the right and left panels 
of Fig. |6] highlights the importance of the absorption from the 
optically thick plasma of the chromosphere in the selection of 
post-shock plasma components that produce observable emis- 
sion in the X-ray band. As a consequence, we found that the 
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X-ray luminosity of the post-shock zone as well as the mass ac- 
cretion rate derived from the emerging X-ray emission strongly 
depend on absorption effects. 

Our exploration of the model parameter space has shown that 
only a fraction of modeled streams predict post-shock plasma 
with measurable X-ray emission. Interestingly, in five out of 
six CTTSs (MP Mus, BP Tau, TW Hydrae, V4046 Sgr and 
RU Lupi) observed with high resolution X-ray spectrographs, 
the pre-shock densities derived from the He-like triplet (10" ^ 
Wacc ^ 10'^ cm"-'), the escape velocities (400 ^ Mace ~ 600 
km s"') and the metal abundances (f < 1) lie in the ranges of val- 
ues that, according to the plots in Fig.|6] are the most suitable for 
the observation of X-ray emission from the post-shock accret- 
ing plasma. The CTTS Hen 3-600 lies at the boundaries of our 
range of observability due to its low escape velocity (Mace ~ 290 
km s"') and the nearly solar metal abundances (^ ~ 1.0). 

Our results can help us to reconcile the (apparent) dis- 
crepancy found betwe en the accretion rates measured from X- 
rays and optical dat a (iSchmittetal .1120051: iGiinther et all 120071: 
lArgiroffi et 311120091: Curran et al. 2010, in prep.). In fact, in the 
case of accretion streams with non-uniform mass density along 
the stream cross section (or in the case of multiple streams with 
different density simultaneously present on the star), we expect 
to observe X-rays preferentially from the low density portion 
of the post-shock plasma, the high density shocked plasma be- 
ing strongly absorbed. For instance, as suggested by 3D MHD 
models of the star-disk system (Romanova et al. 2004), accretion 
streams may consist of a high density central region, surrounded 
by material at lower densities. Consequently, the shocked slab 
should be constituted by a dense core (rooted deeply in the chro- 
mosphere and with a short stand-off height) surrounded by less 
dense shocked material (located shallower in the chromosphere 
and with a large stand-off height). According to our results in 
Sect. [3] the X-ray emission from the dense core is expected to be 
largely absorbed by the stellar chromosphere and by the accre- 
tion stream itself, whereas X-rays emitted from the (less dense) 
shocked plasma at the boundary of the slab may suffer only a 
minor absorption by the chromosphere. All the shocked plasma 
contributes to the optical emission, whereas only a fraction of it 
contributes to the X-ray emission. The result is that the accretion 
rates deduced from optical observations are expected to be larger 
than the rate values deduced from X-ray observations. 

Under the assumption that the accretion streams are inho- 
mogeneous and taking into account the chromospheric absorp- 
tion, our model predicts that different He-like triplets, in general, 
measure different densities of the emitting plasma. In fact, since 
the effects of absorption increases with wavelength, Ne ix (13.45 
A) emission is expected to be less absorbed than Ovii (21.60 
A) emission. Consequently, the high density plasma (i.e. the 
component suffering more absorption effects) should contribute 
more to Ne ix than to O vii emission. In the case of inhomoge- 
neous streams, therefore the average density of the plasma lead- 
ing to Ne IX emission is expected to be larger that that of plasma 
leading to O vn emission. This argument can be checked with 
a simple case of an accretion stream with velocity Mace = 600 
km s"', metal abundance ^ = 0.2 and two components with the 
same cross section, one with density «ace = 10" cm""* and the 
other with Waee - lO'^ cm""*. From Table |3] we obtain that the 
X-ray luminosities in the Ovii resonance line due to the low 
and high density stream components are log(Lovii) = 28.85 
and log(Lovii) = 28.9 erg s"', respectively, whereas the lumi- 
nosities in the Neix resonance line are log(LNeix) = 28.53 and 
log(LNeix) = 29.04 erg s"'. In other words, the denser stream 



contribute for 73% to the Neix emission, but only for 53% to 
the O vn emission. The average density measured with the Ne ix 
triplet is expected therefore to be larger than that measured with 
O vn triplet although Ne ix lines form at temperatures higher than 
O VII lines. 

Finally, our results do not allow us to formulate a unique 
hypotheses to explain why soft X-ray emission from high den- 
sity plasma have not been detected in high mass young accret- 
ing objects (e.g. T Tau, AB Aur and HP 163 296; Gud el et al.1 
l2007HTelleschi et al.ll2007HGunther & Schmitti[2009,) . In partic- 
ular two hypotheses can be formulated: a) accretion streams in 
these objects are, in general, characterized by densities higher 
than Haee = lO'^ cm""*, so that the X-ray emission from shocked 
material is fully absorbed by the chromosphere; b) accretion 
streams are much less dense than streams of low-mass stars 
(«acc < 10'" cm"-'), so that the resulting X-ray emission from 
the shocked material cannot be distinguished from the coronal 
emission. Further analysis of optical and X-ray observations of 
high mass young accreting objects is needed to assess if the lack 
of X-ray emission from high density plasma is a general feature 
of these objects and if some hints on the density of accretion 
streams can be obtained. 



5. Conclusions 

We performed an intensive simulation campaign, using the 1-D 
hydrodynamic model introduced in Paper I, exploring the model 
parameter space which is representative of almost all low-mass 
CTTSs observed to date in the X-ray band. The model describes 
the impact of an accretion stream onto the chromosphere of a 
CTTS under the hypothesis that yS «; 1. Our aims include: 1) 
to investigate the physical properties of shocked slab resulting 
from the stream impact as a function of the density, velocity and 
metal abundance of the accretion stream and 2) to investigate the 
observability of the post-shock plasma in the X-ray band, taking 
into account the absorption from the optically thick plasma of the 
surrounding chromosphere. Our main results can be summarized 
as follows: 

1. The post-shock zone resulting from the impact of the ac- 
cretion stream onto the stellar surface is characterized by 
temperatures ranging from 0.5 to 5 MK (depending on the 
stream velocity) and by stand-off heights of the post-shock 
zone spanning six orders of magnitude from » lO'* to » 10'" 
cm (mainly depending on the stream density). 

2. In all the case the post-shock zone oscillates quasi- 
periodically. Oscillations are composed by a heating phase 
during which the accretion shock builds up a slab of hot 
plasma and a cooling phase during which the post-shock 
zone cools down under the effect of thermal instabilities 
at its base. The oscillation period ranges from w 3 x 10"^ 
to w 4 X lO-' s within the space of physical parameters. 
However, these oscillations are very difficult to observe, be- 
cause the accretion stream is, most likely, inhomogeneous 
and clumped (i.e. with variable accretion rate) and, in the 
case of yS «: 1, constituted by several different fibrils with 
different instability periods and random phases, leading to 
no evident periodic variations of the X-ray emissiorQ 

3. The effect of the absorption from the local chromosphere on 
the X-ray emission from the post-shock zone strongly de- 
pends on the accretion stream properties. The stream density 
is a key parameter to make observable in the X-ray band the 



In the case of /? > 1, see discussion in lOrlando et alj dlOlOh . 
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post-shock material: the higher the stream density, the thin- 
ner and more deep rooted in the chromosphere the shocked 
slab, and the larger the absorption of X-ray emission from 
the slab. We found that, in general, high density accretion 
streams (n^cc > 10'^ cm"^) produce shocked slab which 
are strongly absorbed. Shocked plasma from high velocity 
streams is more easily observable than that from low veloc- 
ity ones, the stand-off heights of the shocked slab increasing 
rapidly with the stream velocity. 
4. Our results suggest that the discrepancy between the mass 
accretion rates derived from optical and X-ray data as well 
as the different densities measured, in general, from the O vii 
and Neix line triplets diagnostics can be explained if ac- 
cretion streams are inhomogeneous (and/or multiple streams 
with different densities are present simultaneously). In fact, 
only the light component of the post-shock plasma is pref- 
erentially observed in the X-ray band, leading to an under- 
estimation of the mass accretion rate with respect to the rate 
value deduced from optical observations. Also, the absorp- 
tion in O VII lines is larger than in Ne ix lines and, therefore, 
the weight of dense plasma is expected to be larger in Ne ix 
than in O vii. 
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